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Abstract 

In this paper a Malliavin calculus for Levy processes based on a family of true derivative 
operators is developed. The starting point is an extension to Levy processes of the pioneering 
paper by Carlen and Pardoux |8] for the Poisson process, and our approach includes also the 
classical Malliavin derivative for Gaussian processes. We obtain a sufficient condition for the 
absolute continuity of functionals of the Levy process. As an application, we analyze the absolute 
continuity of the law of the solution of some stochastic differential equations. 
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1 Introduction 

Over the past twenty five years, it has been done several attempts for extending the Malliavin calculus 
to a scenario driven by processes with jumps. Among these efforts, we should cite the celebrated 
papers of Bismut 0, Leandre [fT8ll and the monograph of Bichteler, Gravereaux and Jacod (5J. In the 
latter it is assumed the restriction that the Levy measure v is absolutely continuous with respect to the 
Lebesgue measure and the main idea is to derive with respect to the size of the jumps. However, if 
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the measure v is discrete, the condition in [0 is not satisfied, and by this reason Carlen and Pardoux 
[8| proposed, for a Poisson process, a derivative operator based on translation of the jump times. 

There is another approach to define a Malliavin derivative as an annihilation operator through a 
chaos expansion. In this way, Nualart and Vives ll23l proved that for the standard Poisson process 
the annihilation operator (Malliavin derivative) agrees with a translation operator. Immediately it ap- 
peared the necessity of extending these methods to Levy processes with a more general Levy measure, 
and many authors proposed and studied such difference operator: see for example, Picard [24] and 
E51 . members of the norwegian school of probabilities like Di Nunno, Meyer-Brandis, 0ksendal and 
Proske [13J and Benth and Lokka JH, and also Ishikawa and Kunita lfT5l . and Sole, Utzet and Vives 
(261. It is a l so interesting to cite the extension to abstract Fock spaces by Applebaum [0. Note that 
in all those works, in contrast with the Wiener case, the annihilation operator is not a true derivative 
and it does not satisfy the chain rule. 

A number of authors have made contributions in the direction of defining a true derivation operator 
for jump processes, under more or less restrictive assumptions. We should cite Denis ifTTTl and Bouleau 
and Denis 0, who use the context of Dirichlet forms, Di Nunno lfl2ll . and recently Kulik [17] and 
Bally Q. We have to note that Decreusefond and Savy IfTOl worked with marked point processes 
(i.e., processes with compensator \(s)dsu(dz)), introducing an operator similar to the one studied in 
this paper, but with the restriction that v should be a probability measure. 

In this paper we define a true derivative operator for functionals of a general Levy processes which 
can have a continuous part. Such operator is an extension to Levy processes of the one introduced by 
Carlen and Pardoux [8 1 for the Poisson process. Our main finding is that this operator can be restricted 
to a Borel subset of E, say A, where live part of the jump sizes of the Levy processes, and a suitable 
(time and space) weight function k(t, x) can be introduced. That subset and weight function can be 
chosen in every application to suit the properties of the problem in hand. By this reason we call it the 
local Malliavin derivative. Our procedure is, as in the Brownian Malliavin calculus, to start with a 
class of smooth random variables, define the derivative, and extend by density to a domain contained 
in L 2 {VL). We observe that this approach can be used even for the difference operator as it is done in 
Geiss and Laukkarinen [14J, but, as we commented, unfortunately, the chain rule is not satisfied for 
the difference operator. 

We would like to point out that we don't present here a complete Malliavin calculus theory, be- 
cause our integration by parts formula (Proposition [377]) does not cover the general case. On the other 
hand, we do not give any result about the smoothness of the density. However we believe that we may 
obtain some interesting results from our development. 

Using the chain rule of our derivative operator we give a sufficient condition for the absolutely 
continuity of functionals of the involved Levy process analogue to the one of Carlen and Pardoux 
[H. In particular, we apply that condition to deduce the absolute continuity of solutions of different 
stochastic differential equations. The first case studied is the equation governed by a Levy process 
with continuous part 

Z t =x+ f b(Z s )ds+ [ a(Z s )dW s + [ [ l(y)ydN(s,y), t e [0,T\, (1.1) 
Jo Jo Jo Jw 

with suitable coefficients b, a and I, and where M := E — {0}. We get a similar result as the one 
contained in Nualart [22J, where jumps are not considered. 

The second application is devoted to the following equation, driven by a pure additive jump pro- 
cess, 

Z t = x+ [ f(Z s )ds+ [ [ h(y)dN(s,y), te[0,T]. (1.2) 
Jo Jo Jm.o 
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Here / : R — >■ R is a function with a bounded derivative and h E L 2 (JH , v) fl L 1 (Mo; z/ ) ? v 

is the Levy measure corresponding to the driving Levy process. It is well-known that this equation 
has a unique square-integrable solution. The first step is to find the derivative of the solution at time 
T. Then, in the case that z/(Ro) < oo and that the function / is monotone, choosing the convenient 
k(t, x), we obtain that Z T satisfies the criterion established for the absolute continuity in the set of 
paths where there is at least one jump. Moreover, if u(Rq) = oo, it is only necessary the condition of 
monotony for / in a neighborhood of the initial condition x. This is the same type of result obtained 
by Nourdin and Simon ll2D using the stratification method due to Davydov et al. flU. Note that as 
the Levy measure is not finite, we have, with probability one, infinity jumps in any neighborhood of 
the origin, and this fact guarantees that we can find the suitable function k in order to restrict us to the 
neighborhood of monotony of /. 
Finally, we study the equation 

Z t = x+ [ f(Z s )ds+ [ [ h(y)g(Z s _)dN(s,y), t E [0,T], (1.3) 
Jo Jo Jm. 

with some conditions on /, g and h. The existence of a density for Z T was analyzed by Carlen and 
Pardoux [8] in the case that the involved Levy process is a Poisson process. When the Levy process 
is a compound Poisson, using our criterion of absolute continuity with a suitable function k(t, x) that 
allows us to concentrate in only one jump, we prove the absolutely continuity of the solution in the 
set of paths where there is at least one jump in the case where a condition on the Wronskian of / and 
g is satisfied. 

The paper is organized as follows. Section 2 presents the preliminaries and notations. Section 3 
is devoted to the study of the family of derivative operators introduced and to the criterion of absolute 
continuity. In Section 4 we present some properties of the derivative operator applied to multiple 
integrals which are needed in the next section. Finally, in Section 5 we get the absolute continuity of 
the solutions to the stochastic differential equations (11.11) . (11.21) and (11.31) . 

2 Preliminaries 

In this section we introduce the framework and the notation that we use in this paper. Let X = 
{X t , t > 0} be a Levy process, that is, X has stationary and independent increments, is continuous 
in probability, is cadlag and X = 0, with triplet (7, a, v), and it is defined on a complete probability 
space (O, J 7 , P); for all these notions we refer to Sato G71l . By convention, the Levy measure v is 
defined on R and ^({0}) = 0. Also, let {F t , t > 0} be the natural filtration associated with X. The 
process X admits the Levy-Ito representation: 

X t = jt + <rW t + / xdN(s,x) + \im / xdN(s,x), (2.1) 

J J (0,t]x{\x\>l} e +° J J (0,t]x{e<\x\<l}) 

where W is a standard Brownian motion, iV is the jump measure of the process defined below, W 
and iV are independent, and dN(t,x) = dN(t,x) — dtv(dx). The limit is a.s., uniform in t in any 
bounded interval. 

For any Borel set B E B((0, 00) x R ), with R := R — {0}, the jump measure is defined by 

N(B) = card{t E (0, 00) : (*, AX t ) E B}, 

where AX t = X t - X t — 
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It is well known (see Ito [ 1 6 1 ) that X determines a centered independent scattered random measure 
M on [0, oo) x R. To this end, consider the measure 

fi(dt,dx) = \(dt)5o(dx) + \(dt)x 2 b , (dx), 

where A is the Lebesgue measure in R. For E E B([0, oo) x R) such that /u(-B) < oo, 

M{E) = a [ dWt+lim / / xdN{t,x), 

JE(0) n->ooJ J {( t ,x)£E':±<\x\<n} 

where £(0) = {* > : (t, 0) G ,5}, E' = E — {(t,0) € £}, and the limit is in L 2 {9). For 
A,Be B([0, oo) x R)), with fi(A) < oo and fi(B) < oo, we have 

E[M(A)M(B)] = /i(AnB). 

Because we work in a finite time interval [0, T], we consider the restriction of M to [0, T] x R. 
For a function fa 6 L 2 ([0, T] x R, //), we can construct the integral (in the L 2 (f2) sense) 

M(h) := / h(t,x)dM{t,x), 

J[0,T]xR 

which is the multiple integral of order 1 defined by Ito lfT6ll . This integral is centered, and for g, h E 

L 2 ([0,T]xR,/i), 

E[M(h)M(g)] = [ hgdfi. 

J[0,T]xK 

We can write this integral as 

M(h) = a[ h(t,0)dW t + [ h(t,x)xdN(t,x). 

JO J[0,T]xR 

In this paper we also need to consider integrals with respect to the non compensated counting 
measure, defined by 



/ h(t,x)x dN(t,x) :— / h(t,x)x dN(t,x) + / / h(t,x)x dtdv(x), (2.2) 

J[0,T]xR </[0,T]xR Jo Jr 

when both integrals in the right hand side make sense. 

Observe that if X has continuous trajectories, then M is nothing more than the independent 
Wiener measure, and if a = and v = 5\ we obtain the standard random Poisson measure. 



3 The derivative operators 
3.1 Carlen-Pardoux derivative 

The starting point of this paper is the following remark due to Leon and Tudor Il20ll about the Carlen- 
Pardoux derivative JSJ: Let iV be a Poisson process of parameter 1, and N t — N t — t be the compen- 
sated Poisson process. Denote by T"i < • • • < T n . . . the jump times of N. Let h E C 1 ([0, T)). Then, 
the random variable 

Y = I h(s) dN s = h ( T n) ~ [ h(s) ds 

JO T n <T ^° 
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is in the domain of the Malliavin derivative in the Carlen-Pardoux sense and the derivative satisfies 
the following equality: 



DtY = j\'{s)(^-l m {s))dN s 



In the next subsection, we extend formally this property to define a Malliavin derivative for a Levy 
process. 

3.2 The local Malliavin derivative 

Henceforth, S denotes the family of all functionals of the form 

f(M(hi), . . . , M(h n )), (3.1) 

where / is in C£°(R n ) (/ and all its partial derivatives are bounded), hi,...,h n G L 2 ([0, T] x R, /j,) 
and, for all x ^ 0, hi{-, x) G C^fO, T]), i G {1, . . . , n}, and d/i* G L 2 ([0, T] x E, /i) where 9 means 
the partial derivative with respect to time. The set S is called the family of smooth random variables, 
and we will prove later that it is dense in L 2 (f2) (see Proposition l3.5l) . 

We will also consider the family /C of all bounded functions k : [0, T] x E Q — >■ K such that they 
and their partial derivative with respect time are in L 2 ([0, T] x R , X x u) f] L 1 ([0,T] x R , A X u). 

Now we introduce the main tool of this paper. Namely, the family of derivative operators for 
functionals of the Levy process X. 

Definition 3.1 Given k G /C, F G S and A G B(M.), we define 

n 

D ^ F = Y^{d l f){M{h 1 ),...,M{h n ))D^ k M{h l ), te[0,T], 

i=l 

where 

D^ k M(h) = l A (0)ah(t,0) + / k(s,y)d a h(s,y)£ ~ h,T\(s))y dN(s,y). (3.2) 

J [0,T] x (AnR ) 1 

We call D^ ,k F the local Malliavin derivative of F. 

Remarks 3.2 1 . The integral in the right-hand side of (13.21) is well-defined. Indeed, in agreement 



with(l2T2 



S - 

k(s,y)d s h(s,y)(- - l [tjT] (s))y dN(s, y) 

[0,T]x(AnR ) 1 



is finite because k is bounded, and 

" T f s 

/ k(s,y)d s h(s,y){- - l [t , T ](s))y dsv{dy) 

J AnRo 1 

is finite by Cauchy-Schwarz inequality. Moreover, D^ k F G L 2 ([0,T] x O, A <g> P). This 
explain why we need to introduce the family /C in the definition of the derivative operator. Also 
in Section [51 we see that this family /C allows to consider applications to stochastic differential 
equations. 
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2. The function k(t,x) is an essential ingredient of the derivative, and it may change from one 
application to another. To shorten the notation, in general we will omit that k in the expression 

D A,k p 

3. ]fh,dh e L 2 {[0,T] x K 0) A <g> v), then h(s, y)/y and d s h(s, y)/y belong to L 2 ([0, T] x E ,/i). 
Hence, 

DfM(h{s,y)/y)= [ k(s,y)d.h{s,v)£ ~ l m (s))dN{s,y), 

</[0,T]x(AnR ) 1 

where the factor y has disappeared in the integral. 

4. Observe that if A = {0}, then D A is the Malliavin derivative with respect to the Brownian part 
of the Levy process X. 

5. If A = {x} for some i ^ with v({x}) ^ 0, we obtain 

D\ x} M{h) = x k{s,x)d s h{s,x){^ - l m {s))dN:, 

where N* is the Poisson process on [0, T] that counts the number of jumps of height x. More- 
over, if the Levy process is the standard Poisson process, x = 1, and if we take k(t, x) indepen- 
dent of the time parameter, we obtain 

D} l} *M(h) = k(l) £ d s h(s, - h,T](s))dNl 

Like we commented in the previous subsection, that is exactly a constant times the derivative 
defined by Carlen and Pardoux in (HI, as it is shown in Leon and Tudor lEOl . 

6. The integration by parts formula allow us to write 



/ / k(s,y)d s h(s,y)y(^--l [tjT ](s))dsv(dy) 

J AnRo Jo 1 



{k(t,y)h(t,y) - - I k{s,y)h(s,y)ds}yu(dy) 



T 

T 



/ / h(s,y)d s k(s,y)(--l[ t ,T\(s)jyu(dy)ds. 

Jo J AnRo 1 



Thus, 



DfM{h) = l A (0)<Th(t,0)+ [ k(s,y)d a h(s,y)y£-l m (8))dN(s,y) 

-/[0,T]xAnR 1 

+ [ {k(t,y)h( y t,y) - ~ f k(s,y)h(s,y)ds}yu(dy) 

J AnRo Jo 

/ h( s >y) d sk{s,y) - l[t,T](s)) yu(dy)ds. 

J AnRo VJ ' 
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7. If ^ A, then D A M(h) is orthogonal to the constant functions. In fact, the Fubini theorem (see 



Proposition [33]) yields 

rT r s 

k(s,y)d s h(s,y)y(— - l [0 , s ](t))dN(s,y)dt 

^[0,T]xA 1 

r T S 

k(s,y)d s h(s,y)y / (- - l M (t))dtdN{s,y) = 0. 

'[0,T]xA JO 1 

The next result is a Fubini type proposition. Its proof is straightforward and omitted. 

Proposition 3.3 Let f e L 2 ([0, T] 2 x R, A® 2 £g> //). TTzen the following integrals are well defined and 

M(f{t,-))dt = M([ f{t,-)dt),a.s. (3.3) 



In the Brownian Malliavin calculus if F is a smooth random variable such that F = 0, it follows 
that DF = 0. However, for Levy processes matters are different. For example, if X is a standard 
Poisson process and we take h — 1, then F = f(M(h)) = f(N T — T), and F = only implies that 
/ is zero on the set N — T. In the next proposition we show the corresponding property. 

Proposition 3.4 Let F 6 S such that F = a.s. Then D A F = 0, A <g> P — a.e. 

Proof. Let F = f(M(hi), . . . , M(h n )). To consider the most general case, assume that G A. By 
definition, 



D A F = J2(dif)(M(hi), . . . , M(hn))<Thi(t, 0) 



— — 

(*) 



n „ 

+ YW)(M(/i!), . . . ,M(h n )) / k(s,y)d s hi(s,y)(^ - l m (s))y dN(s,y) . 

i= l </[0,T]x(AnR ) 1 



— — 



The term (*) is a Malliavin derivative with respect to the Brownian part (see Sole et al. (261). 
From the independence between the Brownian part and the jumps part of X, by fixing the random 
variables corresponding to the jumps part, we have that 



M(hj) = aW(hj) + a 



where W(h) = J Q T h(t, 0) dW t and Cj is a constant; changing the function / we can take C\ = ■ ■ ■ = 
C n = 0. First assume that the Gaussian vector (W(hi), . . . , W(h n )) has a regular covariance matrix. 
Then it has density, and 

f(W(h 1 )...,W(h n )) = 0, a.s. 

implies 

f(xx, ...,x n ) = 0, X n - a.e. 
It follows that all derivatives of / are zero. 
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If the covariance matrix of (W(hi), . . . , W(h n )) is singular, denote by k the rank of that matrix, 
and (reordering) assume that (W(hi), . . . , W(h k )) has a density. Then there are numbers % = 
1, . . . , n, j — k + 1, . . . , n such that 



W(h j ) = Y,t4w(h i ),j = k + l,...,n. 

i=i 

Given the linearity of Ito-Wiener integral, it follows that A-a.s., 

k 

h i = ^2rf h i, j = k + l,...,n. 

i=i 

Consider the vectorial function 

$ = ($!,...,$„) :R k ^R n 

given by 

^ j (x 1 ,...,x k ) =xj, forj = l,...,k, 

and 

k 

^j(x u ...,x k ) = J^/i-Xi, j = k + l,...,n. 

i=i 

Define g := f o G>. Then 

/(W^x) . . . , ^(/i n )) = <?(W^i) . . . , W^)). 
Since . . . , W(h k )) has a density, we have that 

^ ■ i 

By the chain rule, 

Jg JfJ&1 

where J m denotes the Jacobian matrix of a generic function m. Writing h = (hi, ... , h k ) T , we have 

o = Jg h(t) = j f j*h{t) = J2§f. h M 

3=1 j 

and hence (*) is zero. 

The term (**) is treated as follows. Again by the independence between the Brownian and the 
jumps part, we can assume a = 0. For the sake of simplicity we consider the case n — 2. 
Consider the set 

A m = A n {x e R : l/m < \x\ < to}, 

and denote by N Am the restriction of the counting measure of jumps N to the jumps with size in A m . 
That is, for B e B((0, oo) x R ), 

N Am (B) = #{t : (t, AX t ) e B and AX t e A m }, 
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and 

N rn = ^-([CT] x E ) < oo, a.s., 

due to v(h m ) < oo. The jumps of X in [0, T] with height in A m can be ordered, T\ < T 2 < ■ ■ ■ < T N ™. 
Given N™ = k, the vector (Ti, . . . , T k ) has a Dirichlet distribution on [0, T], independent of the height 
of the jumps. The jump size has a distribution given by P{AXT j E A} = v(A)/u(A m ), A E B(A m ). 

In a similar way, define N Ro ~ Am . Since N Am and jV K(, ~ Am are independent, denoting by EA m 
and E K() _ Am the expectation with respect to each jump counting measure, and as F = implies 
F 2 = / 2 (M(/i! ),..., M(K)) = 0, we have that 



fc>0 



= E[F 2 } = E Ro _ Am E Am [F 2 ] = E Mo _ Am | J^nN? = k}E[F 2 \ N? = k] 
Due to the fact that F 2 > 0, for each k, 



f(J2h l (s j ,y j )y 



h x (s,y)ydv(y)ds + M^Iro-aJ : 
h 2 (s, y)ydv(y)ds + M(/i 2 l Ro -A m ) ) ds x ■ ■ ■ ds k du(y 1 ) 



= E Ro _ Ar , 

k 

^h 2 {s jl y j )y j - I I h 2 (s,y)ydv(y)ds + M(h 2 lw -A m )) dsx - ■ -ds k dv(yi) ■ ■ -dv(y k 

where 5)t(A m ) the simplex 

{{t 1 ,x 1 ;...;t k ,x k ) E ([0,T] x A m ) k : t x < ■ ■ ■ < t k }. 

Therefore, using again that f 2 is non negative, P Ro _ Am -a.s., the function from ^(Am) into E, 

k r T r 

{s^yi, . . .;s k ,y k ) ^ f(^2h 1 {s j ,y j )y j - / / h x {s, y)y du(y)ds + M (Mno-A m ) 

j = l JO J Km 



r 



h 2 (s,y)y du(y)ds + M(h 2 l Mo _ Am ) 



is zero a.e. dsi . . . ds k v(dyi) . . . u(dy k ). Hence, by the continuity of both / and ^(-, y), it is deduced 
that, fixed (y%, . . . , y k ), u k -a.e., the derivatives with respect to the time variables s are also zero, and 
then, for every £ = 1, . . . , k 

2 k 



3=1 
k 



J A 

T 



hi{s,y)ydv(y) + M(/iil Ro _ Am ) , 



5^/i 2 (sj,%)% - / / h 2 {s,y)ydiy{y) + M(h 2 l Ro -A m ))d s hi{s i ,y i )y e = 0, 

j= l Jo J A m ' 

a.e. dsi . . . ds k v{dy{) . . . u(dy k ). Multiplying this quantity by k(s£, yg) — l( tjT ](si)) and adding 
with respect to £, we arrive to 



i=i 



[0,T]xA r , 



k(s,y)d s hi(s,y)y(— - l( t) T](s)) dN(s,y) = 0, a.s. 



Now, thanks to Remark 13 .21 we can apply the dominated convergence Theorem to pass to the limit 
that expression when m — > oo, and we obtain that (**) is 0. □ 
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Proposition 3.5 S is dense in L 2 (Vt). 



Proof. Since the finite linear combinations of multiple integrals are dense in L 2 (Vt), it is enough 
to prove that the random variables of the form M(E 1 ) ■ ■ ■ M(E n ) for E x , . . . , E n e B([0,T) x R) 
disjoint sets, with finite ji measure, can be approximated by elements of S. We proceed in four steps. 

1. It suffices to prove the approximation for disjoint closed sets. This is proved in the following 
way. The measure /i is a-finite in <B([0, T] x RJ, and hence it is regular in the sense that for every 
A e B([0,T] x R), 

fi(A) = sup{/i(C) : C C A, C closed} = w£{fi(0) : Ad 0,0 open}. 

So, given Ei,...,E n G B(jO, T] xl), for every e > 0, there are closed sets Ci, . . . , C n , with 
Cj C Ej and 

fi(Ej - < e, j = 1, . . .,n. 
From the fact that Ex, . . . , E n are disjoints and Cj C -Ej, it follows that 



E 



(M(E-l) ■ --M(E n ) - Mid) ■ ■■M{C n )) 2 ] < n max { t i(E j )} n - 1 max {^(E, - CM 

J 1=1, ...,n 7=1, —,n 



and the claim of this step follows. 

2. We can assume that E\ , . . . , E n are disjoint open sets: From previous step, we can take E\ , . . . , E n 
as disjoint closed sets. By the separability properties of W 1 , they can be separated by open sets, that 
is, there are disjoint open sets Oi, . . . , O n such that Ej C Oj. By the regularity of \i, and the same 
reasoning as before, the property is proved. 

3. Let Ei, . . . , E n be disjoint open sets. Since every open set on [0, T] x R is a countable union 
of open rectangles, we can assume that each Ej is a finite union of open rectangles (with finite /i- 
measure) of the form (t',t") x (x',x"), or [0,t) x (x',x") or (t,Tj x (x',x") ; moreover, we can 
reduce slightly the size of all the t-intervals in such a way that the measures f^(Ej), j = 1, . . . , n, do 
not change significantly: specifically, given e > 0, there is S > such that if we denote Ej the union 
with rectangles (£' + 5, t" — 5) ) instead of the corresponding (t',t") x (x',x") of Ej, and 
similarly for the other types of rectangles, then 

fi{Ej-E 5 j)<e, j = l,...,n. 

Let m r E Q°(R), r > a family of mollifiers (see, for example, Adams 01 Page 29]). For a 
rectangle (£' + 5, t" — 5) x (a/, x") write 



fr(t,x) = (m r * l( t '+S,t"-S)j{t) l(x',x"){x), 



where the symbol * means convolution. Collecting all the f r corresponding to the same EL and using 



the fact that the support are included in Ej for r < S, we obtain that a random variable of the form 
M(gi) ■ ■ ■ M(g n ), with g\,...g n with the required properties, and moreover, with disjoint support, 
approximates the product M(Ei) ■ ■ ■ M(E n ) in L 2 (Q). 

4. Finally, use an approximation by Q°(IR n ) functions to the product function to get the Proposition. 

□ 



Remark 3.6 In the preceding proof we have proved that instead of S we can use the sum of the 
smooth random variables of the form f(M(hi), . . . , M{h n )) with h\, . . . ,h n with disjoint support to 
define the local Malliavin derivative. 
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Next proposition gives a Malliavin integration by parts formula: 
Proposition 3.7 Let F 6 S, and g be a measurable and bounded function on [0, T]. Then 



E 



T 



{D?F)g(t)dt 



E 



1 

T 



g{t)dt)k(s } y)dN(s,y) 



F[l A (0)a I g(s)dW s + / (g(s) - 

</[0,T]x(AnR ) 

2 F( I d s k(s,y)[ [ T g(t)£-l [0>B] (t))dt]dN(s,y) 

[0,T]x(AnR ) JO 1 



(3.4) 
(3.5) 
(3.6) 



Proof: 

First we consider the case F = f(M(h)), where / and h are as in (13.11) and assume v(A) < oo. 
Then, 



E 



(D?f(M(h)))g(t)dt 



1a(0)E 

+E 



a / f'(M(h))h(t,0)g(t)dt 



f(M(h)) 



'0 V[0,T]x(AnR ) 

From the integration by parts formula in the Wiener case, 



k(s,y)d 8 h(s,y)y(— - l m (s))dN(s,y) g{t)dt 



(3.7) 



E 



D$ 0} f(M(h))g(t)dt 



Uo 



E[f(M(h))M(g(t)®l {0} (x))] 



and the first term of the right hand side of (13.71) is equal to 

l A (0)E[f(M(h))M(g®l An{0} )}. 
Now, using Fubini theorem, the second term can be written as 



E 



f(M(h))[ I {k(s,y)d s h(s,y)y(^ / g(t)dt - / g(t)dt)\dN(s,y) 

[0,T]x(AnR ) L 1 JO Jo J 



which, by the independence of W, N A and N R ° A (using the same notations as in the previous proof), 
and the symmetry of the function Y^j=i h{ s ji yj)Vj> ls equal to 

oo 1 



n=l 



u(A) n T r ' 



x / / 




J[0,T]" J(j 


vnR ) n ^ 



h(r, z)z drdu(z) 



IfCEK^y^- [ [ 

1 ~[ Jo Jk 

+ l A (0)M(/il {0} ) + M(/ l l Ro _ A ))] 

Y^H s jiyj) d M s jiyj)yjy7f; J g(t) dt - J g{t)dtjjds 1 ---ds n dv(y 1 )---di>(y n ) 
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Integrating by parts each integral 



T 



E 



W,N R o~ A 



f'C^2Hsj,yj)yj 



3=1 



J A 



h(r, z)z drdv{z) + ■ 



x k(si,yi)d s h(si,yi)yi 



T 



T 



g(t)dt 



g{t)dt\ds 



and reordering terms, we arrive to 

E 



f(M(h)) I k(s,y)[g(s) 

[0,T]x(AnR ) 



T 



g(t)dt)dN(s,y) 



E 



f(M(h)) 



d s k(s,y)( I g(t)(--l [0iS] (t))dt)dN(s,y) 



'[0,T]x(AnR ) 

which implies that the result holds in the case that z/(A) < oo. 

To deal with the case z/(A) = oo, let A n = {0} U < \x\ < n}. Then, changing A by A n fl A, 
and letting n — > oo, we obtain the result. 

Now, consider F = f(M(hi), . . . , M(h n )), with h±, . . . ,h n with disjoint support. Since the 
random variables M(hi), . . . , M{h n ) are independent, we can compute the expectations in (13.41) iter- 
atively, and the result follows. Finally, by linearity, the equality (13.41) is extended to a sum of smooth 
random variables of the form f(M{hi), . . . , M{h n )), with hi, . . . , h n with disjoint support, which 

□ 



generate S, in agreement with Remark [3761 



The following result is a consequence of Proposition l3.7l applied to FG and of the fact that D A is 
a true derivative operator, which follows from Definition l3.ll 

Corollary 3.8 Let g measurable and bounded, and F,G G S. Then 



2 ( C l T {DtF)g{t)dt + F f T (DfG)g(t)dt 



'o 

E 



FG(\ A (0)^ ag(s)dW s + ^ 



[0,T]x(AnR ) 



g{s)-- I g{t)dt)k{s,y)dN{s,y) 



d s k(s,y) 

[0,T]x(AnR ) V<i 



s \ 
g(t)(--l [0>s] (t))dt) dN(s,y) 



Now we can proceed as in Nualart Il22ll (pag. 26) to show the following result, taking into account 
that the bounded functions are dense in L 2 ([0, T]). 

Corollary 3.9 The operator D A ,from L 2 (Q) into L 2 ([0, T] x Q), is unbounded, densely defined and 
closable. 

In particular, we have that the operator D A has a closed extension, which is also written by D A . 
The domain of this operator is denoted by DomD A . Note that F belongs to DomZ) A if and only if 
there is a sequence {F n : n > 1} C S such that F n converges to F in L 2 (f2), when n — > oo, and 
{D A F n : n > 1} converges to some process Y in £ 2 ([0, T] x fl). In this case, we have that_D A F = Y. 

The following rules of derivation are well known in Brownian Malliavin calculus; they can be 
translated to our context with the same proof. 
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Proposition 3.10 1. (Chain rule). Let f : M. n ~&L be a continuously differentiable function with 
bounded partial derivatives, and let F = (Fx, . . . , F n ) a random vector such that Fj G Dom D A , for 
j — 1, . . . ,n. Then f(F) G DomD A and 

n 

D A {f(F)) = Y,d 3 f(F)D A F r 

3=1 

2. Let F,G e Dom D A such that G D A F and F D A G belong to L 2 ([0,T]xQ). Then F G G DomD A 
and 

D A (FG) = GD A F + F D A G. 
3.3 Absolute continuity via the operators D A 

The chain rule allows to prove a criterion for the absolutely continuity of a functional F in a similar 
way that in the Brownian case. Specifically, the result obtained by Carlen and Pardoux [8] Theorem 
4.1] for the Poisson process can be extended to the Levy case; the proof is the same as Carlen and 
Pardoux, and omitted. 

Proposition 3.11 Let A G B(R), k G /C and F G DomD A such that 

f [DfFfdt > 
Jo 

a.s. on a measurable set A 6 J 7 . Then, P o F^ 1 is absolutely continuous on A, that is, X(B) = 
implies P({F G B} n A) = 0. 

It is worth to remark that the criterion is true for every set A and weight function k G /C. This is very 
interesting for applications because we can choose an appropriate A and k depending of F. It will be 
used in Section \5\ to see that the solutions of some stochastic differential equations have a density on 
a measurable set. 

4 Multiple integrals over sets with finite Levy measure and its 
derivatives 

In this section we study multiple integrals defined on a set with finite Levy measure. Such integrals 
can be computed pathwise, and its properties can be proved by combinatorial methods. 

Let G B(R) be a bounded set such that 0^9 (closure of 0); then u(Q) < oo. As we have 
done before, consider 

N e (B) = #{t : (t, AX t ) G B and AX t G 0}, for B G B((0, oo) x R ), 

that is a Poisson random measure with intensity A <g> vq, where vq(C) = v(Q fl C), for C G <6(R). 
Write 

N® = N e ([0,T] x 0) < oo, a.s. 
Then we can order the jumps in the interval [0, T], say T\ < • ■ ■ < T N &, and write 

3=1 
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Let <p : [0, T] x ]R — )• IR be a measurable function (indeed, it is only needed that <fi is defined in 
[0, T] x 0) integrable with respect to N e and consider 

Jf(0) : = I 4>(t,x)dN (t,x). 

J[0,T]xR 

This integral has a compound Poisson distribution with characteristic function (see, for example, Sato 
11271 Proposition 19.5]) 

E [ e ^i e W] = eX p { f _ l) dtdu(x)}. 

J[o,T]xe 

Thus, if G L p ([0,T] xl ,A® i/ e ), then J®(0) G L p (tt), forp > 1. 

Also multiple integrals can be considered. Recall the notations introduced in the proof of Propo- 
sition [331 let S n (Q) be the simplex 

{(t 1 ,x 1 ;...;t n ,x n ) G ([0,T] x9)": < h < ■ ■ ■ < t n < T}. 

For0 : S n (Q) ->• R, define 



0(ti, xi; . . . ; t n , z n ) dN ^, x x ) ■ ■ ■ dN {t 

m ■En) 



tBi 



S n (B) 



cl ) (T h ^X Th ;...;T in ,AX T Jl 



{TVf >n}- 



l<ii <■■■<«„ <N® 

Proposition 4.1 For every p > 1, i/0 G L p (S n (6), (A ® z/ ) n ), ^en J®(0) G L p (fi), and 

E[|J„ e (0)| P ] <C p>n [ \cf>(t 1 ,x 1 ;...;t n ,x n )\ p dt 1 ---dt n du{x l )---dv{x n ), (4.1) 

■/Sn(6) 

where the constant C n , p depends only on n and p. 

Proof. The main tool of the proof is that the jump times of a compound Poisson process follow a 
Dirichlet distribution, independent of the size of jumps, which are also independent and each of one 
has law P{AX Tj G A} = u(A)/u(Q), for A G B(Q). To simplify the notations, we assume T = 1. 
By Jensen inequality, 

\Jn(<f>)\ P < ( (V (6) ®- ■ •®iV e )(5 n (6)))^ 1 / Wh, Xl ;...;t n , x n ) | p dN e (t u x t ) ■ ■ ■ dN e (t n , x n ). 
v J Js n (e) 

Hence, 



k>n 



-1/(0) 



E 

k>n 

E 

k>n 



fc! 



i/(6) n 



E |0(T Jl) AX Tn ;...;T ln ,AX Tin ; 



K£i<— <in<fc 



fc! 



u fQ\n E 

v 7 l<ii<— <t 



,<k 



E 



9" 



|0(C/ (il ),a;i; . . . ; U( in ), x n ) \ P du(x 1 ) . ..dv(x n ) 
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where (Uu^, . . . , Uu n \) is the ordered . . . , statistic from a sample U\, . . . , of a uniform 
law in (0, 1). Extend to the whole (0, l) n by symmetry, and denote by C\ the family of all subsets 
of {1, ... , k} of size n. Fixed x\, . . . , x n E 0, 



E 



ii)) • • • i ^(i„)> x n ) 

l<il<---<«n<fc 

n 



E 



) •En) 

{h,...,i n }ec n 



<l)(ti,xi,. . . ,t n ,x n )\ p dti ■■■dt n 



(0,1)" 



k\ 



(k-n)\ J 0<tl ... <tn 



\4>(t 1 ,x 1 , . . . ,t n ,x n )\ P dt ± ■ ■■dt r . 



Finally, 



< / / \(j)(t 1 ,x lj ... J t nj x n )\ P dt l ---dt n du(x l ) ...du(x n ) 

Jo<ti-<t„<i Je™ 



x e 



-,(8) "(e)*- w 

{k-n)\ 



E 

k>n 



kP 



-1 



and the result follows. □ 

Corollary 4.2 Let (fi k , <p e L p (S n (Q), (A ® z/ e ) n ), smc/i ?/za? 

Iim0 fe (ti,a;i; . . .;t ) = Xi; . . . ; t n , i n ), uniformly. 

k 

Then 

lim J n (0 fc ) = ^n(0), 

A crucial result in multiple integration theory is the product formula (see Lee and Shih [ 19 1 and 
Tudor ESI ). For this particular case, the formula is simpler and the proof straightforward: 



Proposition 4.3 Consider 4> n : S n (Q) — > R and <f>\ : [0, T] x R — > R- Tfoen, 



J n e (0 n )^i (0i) = ^ 



(4.2) 



where 



n+l 



<t>n®<t>i{h, a?i; • • • ; Wii ^n+i) = E ^(^1,^1; • • • 

where tj , x j means that this pair is missing, and 

n 

, Xi, . . . , t n , X n )(f>i(tj, Xj). 
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Remark 4.4 Assume that is defined in all ([0, T] x 0) n . The function 0„®0i is not the standard 
symmetrization of <f> n <g> 0i, which is given by 

where the sum ranges over the set of permutations of n + 1 elements. We have the relation 

4>n <8> 0i = 4>n®4>ii where n is the symmetrization of 6 n . 

n + 1 

Lemma 4.5 Let G <B(R) £>e a bounded set such that 0^6, and A G <6(R) wzY/i 6 C A. Consider 
f G C£°(R). 77w?n f(N$) e DomD A and 

D A f(N°) = 0. 

Proof. Write 

iV|= / dA e (t,x) = / -l e (x)xdA(t,x) 

i[0,T]xK i[0,T]xMo ^ 

1 ,T 



— le(x)x ciA(£, x) + / / dtdu{x). 

[0,T]xR X Jo J® 

Hence, A® G 5, and D A A® = 0. Thus also /(A®) G 5 and D A f(N®) = 0. □ 

Given a set A C N, the random variable l^ N e eA y can be conveniently approximated by functional 
/(A®) G Dom Z} A . From the closability of D A we get the following corollary: 

Corollary 4.6 For every set AcN, l{jvf eA} £ DomD A and D A l| Ar e gj4 } = 0. 

Theorem 4.7 Let and A as z>z Lemma 14.51 and G L 2 (S n (Q), (A £g> z^) n ) smc/i that for every 
(x\, . . . , x n ) G n , 0(-, Xi; . . . ; •, x„) G C°° (S n ) , where S n = {0 < ti < ■ • • < t n < T}. Then 
J®(0) G DomD A and 

n 

A A 4 e (0) = E ^(Ksj^ds^xx; . . . ; - l (t ,T](%))). (4.3) 

Proo/ We proceed by induction. For n — 1, as in the previous lemma, write 

J A (0) = / 4>(s,x)dN e (s,x) 

</ [0,T]xR 

— M(1q(x)(P(s,x)/x) + / / 0(s, x) dtdis(x). 



Jo Je 

Its derivative is 

s 

l e (x)k(s,x)d s 4>(s,x)[— - l(t,T](s)) dN(s,x) 

[0,T]xR 1 



which coincides with (14.31) . 
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Now, consider n G (S n (Q)) and X G C£° ((0, T) x i? ) . Then the result is obtained using 
the product formula (14.21) 

4 6 +l(0n®0l) = 4 e (0„) Jf(0l) " J n 9 (0n * 0l), (4.4) 

checking that all terms in the right hand side are in DomD A and using the induction hypothesis. 

For a general G L 2 (S n+ i(Q), (A Cg) z/) n+1 ) with the conditions of the theorem, first extend this 
function by symmetry to ([0, T] x 0) n Since is a bounded set by hypothesis, using a mollifier 
family m £ G Cg°(M n+1 ), define 



4> £ {ti, x\\ ■ ■ ■ ; t n+1 , x n+1 ) = / yi ; . . . ; t n+ i, y n+ i)m £ (a;i - j/i, . . . , x n+ i - y n+ i) 

</]R' l + 1 

■ ■■dp(y n+ i). 

Then, e G C^°(M 2n+2 ), lim £ ^ 0e = 0, uniformly and lim e ^ d tj 4> £ = d tj 4>, uniformly. Hence, by 
Corollary 14.21 lim e ^ J n+1 (0 e ) = J„, + i(0) in L 2 (fi), and also the right hand side of (14.31) correspond- 
ing to e converges in L 2 ([0, T] x fi) to the corresponding to 0. Thus, it suffices to prove the formula 
(14.31) for G Q°(IR 2ri+2 ). For such we can apply Theorem 1, page 65 of Yosida [|29ll , that says 
that we can approximate by finite sums Y%=i V?^ where 4^ e C£°(IR 2n ) and 0j i) G C^(R 2 ), 
and also the derivatives of the sums converge uniformly to the derivatives of 0. Indeed, that theo- 
rem by Yosida E9l shows the result for complex valued functions, however, it can be applied to real 
values functions and, in that case, the approximating functions are also real. From Corollary 14.21 
the symmetrization property of Remark 14.41 and the previous induction step, the proof is complete. 
□ 

Corollary 4.8 Let O, A and be as in Theorem \4.7\ With the previous notations, 

0(T 1; AX Tl ; . . ,r ni AI r Jl {Jv| > n} G DomD A 

and, over {N® > n}, 

A A 0(T 1; AX Tl ; ...;T n , AX Tn ) = £ k(Tj, AX Tj )d^(T u AX Tl] . . . ,T n , AX Tn ) - l [0 , Tj] (t) 

3=1 

Proof: We only need to observe that 

oo 

{N£>n} {N£=£}, 

where 

i>t(si,yi; ■ ■ ■ ; s t , ye) = 0(si, yr,...] s n , y n ), £> n, 



and apply Corollary 14.61 and Theorem 14 .71 □ 

Remark. Note that this result yields that with k — 1, agrees with the operator introduced in 
Carlen and Pardoux [8] for the Poisson process. Consequently, this paper extends the approach given 
in (H to Levy processes. Also, in general, with the previous hypothesis, the jump time Tj is in the 
domain of D A and 

A A T, = k(Tj, AX Tj ) - l m] (t)j . 
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5 Absolutely continuity for stochastic differential equations 



In this section we use Proposition 13.1 II and Theorem 14.71 to find conditions that guarantee that the 
solution at time T of some stochastic differential equations driven by a Levy process has density with 
respect to the Lebesgue measure. The key point is that we can choose the convenient set A and weight 
k(t, x) for each type of equation. 

5.1 An equation driven by a Levy process with continuous part 

Here we consider the solution of the stochastic differential equation with an additive jump noise and 
a Wiener stochastic integral of the form 

Zt = x + [ b{Z s )ds+ [ a(Z s )dW s + [ l{y)ydN{s,y), te[0,T], 

JO JO i(0,i]xMo 

where I E /C and the coefficients b and a are differentiable on R with bounded derivatives. Then, 
choosing A = {0} and k E JC, and proceeding as in Nualart [|22l (Section 2.3), we can see that 

D\ Q} Z T = a(Z t )exp (j\'{Z s )dW s + j* |fe'(Z s ) - ^'(Z,)) 2 } dsY t E [0,T]. 
Thus we can state the following theorem: 

Theorem 5.1 The random variable Z? is absolutely continuous with respect to the Lebesgue measure 
on the set {S < T}, where 

S = inf it E [0, T] : J l {a{Zs) ^ } ds > J A T. 

Remark As in ll2~2ll . we can see that last theorem also holds in the case that the coefficients are only 
Lipschitz functions with linear growth. 

5.2 A pure discontinuous equation with a monotone drift 

In this section we consider the following equation with an additive jump noise 

Z t = x+ [ f(Z 8 )ds+ [ h(y)dN(s,y), tE[0,T], (5.1) 

Jo J[0,t]xR 

where a;6K, / : R — >IRisa function with a bounded derivative and h E -^ 2 (IR , v) fl L 1 (IR , v). It 
is well-known that equation (15.11) has a unique square-integrable solution. 

In order to calculate D R °Z T , we need the following result. 

Proposition 5.2 Let Q m = {x e R : 1/m < \x\ < m}, and 

Zl m) =x+ [ f(zW)ds+ [ h{y)dN{s,y), te[0,T]. (5.2) 

Jo J[0,t]x@ m 

Then, as m — > oo, zf"^ converges to Z t in L 2 (Q), for every t E [0, T\. That convergence is also a.s. 
for every t E [0, T] a.e. 
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Proof: Let M > be a bound for the derivative of the function /. Then 



E 



zl m) - Z t 



^ < 2M 2 T jf E (JZ s (m) - Z a | 2 ) ds + 2E M jj" 



h(y)dN(s,y) 



[0,t]x(Ro-e„ 



The result is a consequence of Gron wall's lemma and of the fact that 



E 



[0,i]x(R -e m ) 

E 



h{y)dN{s,y) 

h(y)dN(s,y) 



[o,t]x(R -e m ) 



< 2 



h(y) 2 u(dy)ds + 2 



./Ro-©* 
T 







h(y)v(dy)ds 



\h(y)W(dy)ds 



goes to zero as m — > oo because 6 m j" K - Finally the second claim of this result follows using a 
similar argument. Thus, the proof is complete. □ 



Now we follow a similar reasoning as Carlen and Pardoux (8]|. Consider the flow {$t(s, x) : < 
s < t <T and igR} associated with equation (15.11) . It means, $t(s, a;) is the unique solution to the 
equation 

$i(s,x)=x+ / f($ u (s 7 x))du, te[s,T\. 

J s 

Note that the partial derivatives of the flow with respect to time is 

d t ^t(s, x) = f($ t (s,x)). 
The partial derivative with respect the initial condition satisfies the equation 



d x $ t (s,x) = l+ I f'(<S> u {s,x))d x <5> u {s,x)du 

J s 



and so 



d x <5> t {s,x) = exp QT f($ u (s,x))diij . 
Similarly we have 

a $ t (s,z) = -/($ s (s,z))+ / / / ($„(s,x))9 s $ tt (s,x)du. 

J s 



Hence 



d s <5> t {s,x) = -/(x)exp ^ /'($ tt (s,a;))du^ . 



Denote by X = {X t , t > 0} the Levy process associated to the Poisson measure N as in (|2.1I) 
(obviously with 7 = a = 0). The processes X and Z jump at the same times, and the jumps height of 
the solution process Z in a jump time r is h(AX T ). Denote by Ti < T 2 < .. the jump times of N Qm 
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(the dependence of Tj on m is suppressed to short the notations). The solution of (15.21) in {N® m = n} 
is given by 



47 ] =$ Tl {0,x) + h{AX Tl ), 



and so on. See Figured! 




Figure 1. A trajectory of Z t 



(m) 



By Corollary 14.81 and induction over n 



•)R y(m) 



f(Z ( Ak(T n ,AX Tn 



T 



1 [0,T„](*) 



i=l 



f(4 m) H (7(4™-) - w,Ax r j (§-i [0)Ti] (t) 



on the set |iV® m = n}. Thus we have that 



n R / y(m) 1 



l{JV® m =n} 



/(4: ) )exp Qf /'(ZW)d«j fc(T B , AX T J - 1 [0) T„](0 



,(m)- 



-{Af® m =n} 



£expQf /'(4 m) )^) (/(4^)-/(4 

f{Z^)du 

x (/(4- } ) - /(^ Am )) k(8,y) (~ - l [M (t)) dAT( 8l y). 



x/ 



1 {N^=n} I eX P 

'(0,T]x9 m \./.s 



(5.3) 
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Now we can see that the random variable Z T is in the domain of D R °' k . 
Proposition 5.3 Under the assumption of this section, we have that £ Dom D R °. Moreover, 

Df°Z T = [ exp ( f f'(Z u )du) (f(Z s .) - f(Z s )) k(s,y) - l M {tj) dN(s,y). (5.4) 

J[0,T\xR \Js / \1 / 

Proof: The result follows from Proposition l5.2[ letting m — > oo and from the dominated convergence 
theorem. That theorem can be applied due to the integrand in the right-hand side of (15.31) is bounded 
by M\h(y)k(s, y)\, for some M > 0, since /' is bounded, and the jumps of the solution are related 
with h. □ 



5.2.1 Case with finite Levy measure 

In this section we analyze the existence of a density for Z T in the case that the Levy measure is finite 
and the drift is a monotone function. 

Theorem 5.4 Assume that z^(IRo) < oo. Consider equation 0.71) . where f is a momotone function 
with bounded derivative and h £ L 2 (M. ,u)) with h(y) ^ for y £ R - Then Z? is absolutely 
continuous on the set {/V* > l}. 

Proof: We first assume that / is an increasing function. We choose for k the function (— h A 1) V (— 1) 
that depends only on y. Then, (/(Z s _) — f(Z s ^ + h(AX s )))k(s, AX S ) is strictly positive in the jump 
points. Thus, to finish the proof, we only need to apply Proposition 13.1 H and observe that 

cT 



/ 2 
(|;-W*)) dt = s(l-^), se[0,T], 



and 

rT 



/ (| - l[o, a] (*)) (f ~ l[o,r](*)) dt = s(l - £), forO < s < r < T. 

J 

Hence the two integrals are positive for s < t. Indeed, we have that 

r (Df°z T ) 2 dt 

Jo 

N T° ( T \ 

= E ex P f 2 J Tf0 f'( Z u)du\ (f(Z Tf0 -) - f{Z Tf0 )f (k(AX Tf0 )) 2 T^ (l - 
+2 ex P ( jC f'( Z u)du) (f(Z Tf0 -) - f{Z Tf0 )) k{AX^ ) 

jT f{Z u )d^\ (/(^o-) - f{Z^ )) k{AX T ^ )Tf» ^1 



T 



l<i<j<N T 

x exp 



rpK 

3 



which is bigger that zero on the set {/V* > l}. Therefore Zt absolutely continuous in this set. 

Finally we can proceed similarly in the case that / is decreasing using the function k = (h A 1) V 
f-D. □ 
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5.2.2 Case with infinite Levy measure 

Now we deal with the case that / is only monotone on an neighborhood of the initial condition x. 
This problem has been analyzed in Nourdin and Simon ETTl using a stratification method. 

Theorem 5.5 Consider equation d5.il) . where f is a function with a bounded derivative and h G 
L 2 (IRo ) z/ ) H L 1 (]Ro ; v)- Assume that i/(Ro) = oo, h(y) ^ Ofor y G Mo. and that f is monotone on a 
neighborhood of the point x. Then, the random variable is absolutely continuous. 

Proof: Here it is convenient to write the function k that we use in the derivative, D A ' k . Proceeding as 
in Theorem 15 .41 we assume that there is e > such that / is increasing on (x — e, x + e); the proof 
for / decreasing is similar. 

Note that (15.11) . together with the fact that / has a bounded derivative, implies that there exists 
M > such that 

\Z t -x\<M f \Z s -x\ds + t\f(x)\+ f \h(y)\dN(s,y), te[0,T\. 

JO J[0,t]xR 

Consequently, Gronwall's lemma leads to write 



\Z t -x\<(t\f(x)\+ \h(y)\dN(s,y))e MT , te[0,T\. 

V J[0,t]xM. o J 

Therefore, there is t G (0, T) such that 

\Z t -x\< £ - + e MT [ \h(y)\dN(s,y), te[0,t }- (5-5) 



Now, fort G [0,t ], let 



A t ={e MT f [ \h(y)\dN(s,y)> £ - 
Jo Jr a 



and kt ■ [0, T] — > IR a function in C 1 ^^)) such that k t (s) > for s G [0,t), and k t (s) = for 
s > t. Define k®(s,y) = —k t (s)((h(y) A 1) V (—1)). As the Levy measure is oo, there are, with 
probability 1, jumps in (0, t), and from (15.51) we deduce that the jump sizes satisfy \Z S — x\ < e for 
s < t < to on the set A\. 

Indeed, in this case we have that, on A c t , {f{Z s J) - f(Z s ))k®(s, y) > 0, for < s < T, with 
(f(Z s _) - f(Z 8 ))k®(s,y) > if Z has a jump at times G (0,t). So. (1541) gives 

£ (D*o,kU ZT y du > o in A c t . 

As a consequence, by Proposition 13.1 1[ the condition X(B) = implies that 

P({Z T G B} n A c t ) = 0. 
Finally, we observe that A c t C A\, for t' < t, and using the Markov inequality, we have 

2MT ( / r \ 2N 



P(A t )<^^E[( [ \h(y)\dN(s,y)) 



8e 2 



[0,t] xM 

2MT ( r / r \ 2 

|2. 



22 



Hence we see that lim^o P(A t ) = 0. 

So P(U t<t(ht€Q A c t ) = 1, and we can conclude, as P{{Z T G B}) = \mv m P({Z T G B} n A c t ), 
that Z T is absolutely continuous. □ 

5.3 A pure discontinuous equation with non zero Wronskian 

In this section we assume that the Levy measure v is finite and consider the stochastic differential 
equation 

x + f f(Z s )ds+ f [ h(y)g(Z s _)dN(s,y), t G [0,T], (5.6) 
Jo Jo Jr 

where x G R, h ^ is a bounded function in L 2 (IRo? z/ )» an d / : R — > R and g : R — >■ R are 
two bounded functions with two bounded derivatives and one bounded derivative, respectively. The 
existence of a density for Zt was analyzed by Carlen and Pardoux §8§ in the case that the involved 
Levy process is a Poisson process. In the remaining of this section we also assume that 

\h(y)W(g,f)(x)\ > ^llf HoolWlLlMlL x G R y G R , (5.7) 

where W(g, f) = g'f — f'g is the Wronskian of g and /. 

As in Section 15^21 we consider the flow {<&t(s, x) : < s < t < T and x G M} associated with 
equation (15.61) and the family of jump times {Tj : j G N} of the Levy process X. Remember that 
now we are dealing with a finite Levy measure in this part of the paper. 

Note that the fact that Z Tn = $ Tn (T n _i, Z Tn _ x ) + h(AX Tn )g(Z Tn _) jointly with T = and 
Corollary 14.81 allow to utilize induction on n to see that, for any k G K, 

Df°Z Tn = f(Z Tn ^k(T n ,AX Tn )f^-l [0jTn] (t) 

+f(Z Tn _)h(AX Tn )g'(Z Tn „)k(T n ,AX Tn ) - l [0 ,T n ](*)) 

Vexp ( f n f{Z u )dvt)r t {T u AX Ti ) 



fT, 



+ 5Z / exp ( / f'( Z u)du\h(y m )g'(Z Sm _) ■ ■ ■ h{y 2 )g'(Z S2 _)r t (s 1 ,y 1 ) 

m=2 J Srn(T m ) \J s 1 J 

x dN(s m , y m ) ■ ■ ■ dN(si,yi), (5.8) 



where 



and 



r t (s,y) = k(s,y) - l M (t)) (f{Z a -) - f(Z s ) + f{Z s J)h{y)g'{Z s J)) 



S m (r) = {( Sl , yi ; s m ,y m ) G ([0,T] x R ) m : < s, < ■ ■ ■ < s rn < r}. 
Now we are ready to see that Z T is in the domain of D R ° . 

Proposition 5.6 Let Z be the solution of equation (\5.6\) . Then Z T G DomD R °. 
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Proof: By Corollary 14. 8 1 and (15.81) . we have Z T l^ N n _ n ^ is in the domain of D and 

Df°(Z T l [N ^ =n] ) = Df°($ T (T n ,Z Tn )l {<0=n} 

= 1 [iV*0=r 



■/(Zrjexp /'(Z u )dri)*(T B) AX r J - l [0 ,T n ](t)) 
+ exp ^ f'(Z u )du^D*>(Z Tn 



1 [N^=n] \zZ exp ( f'(Z u )du\r t (T h AX Ti ) 

i<n \JTi J 

+ Yl / exp ( / f'( Z ^) du ) Kym)g\Z Sm -) ■ ■ ■ h(y 2 )g'{Z S2 ^)r t {s 1 , y x ) 

x dN(s u yi) ■ ■ ■ dN(s m , y m ) J . 
Hence, the fact that there is a constant C > such that \r t (s, y) | < Ck(s, y) and 
/ exp I f'(Z u )du ) 

JS m (T n ) \J Sl J 

x h{Vm)g'{Z Sm J) ■ ■ ■ h(y 2 )g'(Z S2 -)r t (s 1 ,y 1 )dN(s 1 ,y 1 ) ■ ■ ■ dN(s m , y m ) 

(HP \ ✓ \ Tit — 1 

f [ \k(s,y)\dN(s,y)) ( f [ \h(y)\dN(s,y)) 
JO JM.Q J \Jo JR J 

implies that Z T is in the domain of D R ° with 

Df°Z T = [ expff f'(Z u )du)r t {s,y)dN{s,y) 

J[Q,T]xR \Js J 

+ Y1 / exp ( / f i y Z u)du)h(y m )g\Z Sm ^)---h(y 2 )g'{Z S2 _)r t (s u y l ) 
m=2 Js m (T) \J Sl J 

xdN(s 1 ,y 1 )---dN(s m ,y m ). (5.9) 

Thus the proof is complete. □ 

Finally we can see that Zt has a density. 

Theorem 5.7 Let Z be the solution of equation ( 15.(51) . Then Zt has a density on the set [N^° > 0]. 
Proof: For a rational number p E (0, T) and a positive integer n > 1, we introduce the set 

A p , n = {N T = n}n {T„_x < p < T n }, 

with the convention T = 0. We choose a function h p : [0, T] — > R+ of class C 1 ((0,T)) such that 
/i p (s) = for s < p, and h p (s) > for s > p. Then, (15.91 ) implies that in the set A p>n we only have 
to consider in the derivative the jump at T n , 

Df» h "Z T = exp f\Z u )d u y p (T n ) - l [0 , Tn] (t)^ 

x (/(Zr B -) - /(Z T J + f(Z Tn .)h(AX Tn )g'(Z Tn .)) , 
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which is different than zero due to condition (15.71) . Hence Proposition 13. 1 1 1 implies that P({Z T E 
B} n A p>n ) = for any Borel set B e B(R) such that X(B) = 0. Thus, for this Borel set we have 

P({Z T E B} fl {N T > 0}) = P({Z T £ B} H (U neN , peQ A p , n )) = 0, 

and the proof is finished. □ 
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